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Quantum dots with large Thouless number g embody a regime where both disorder and interac- 
tions can be treated nonperturbatively using large- A'^ techniques (with N = g) and quantum phase 
transitions can be studied. Here we focus on dots where the noninteracting Hamiltonian is drawn 
from a crossover ensemble between two symmetry classes, where the crossover parameter introduces 
a new, tunable energy scale independent of and much smaller than the Thouless energy. We show 
that the quantum critical regime, dominated by collective critical fluctuations, can be accessed at 
the new energy scale. The nonperturbative physics of this regime can only be described by the 
large- A'^ approach, as we illustrate with two experimentally relevant examples. 
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PACS numbers: 73.50.Jt 

Mesoscopic systems show many novel properties 
which are not present in bulk systems [ij, such as the 
Coulomb Blockade in the zero-bias conductance through 
a quanutm dot, or persistent currents in mesoscopic rings 
penetrated by a magnetic flux. Disorder renders the 
single-particle states chaotic [2 , with the statistics of the 
energies and wavefunctions controlled by Random Matrix 
Theory (RMT)|3|, which also controls the correlations 
between different states separated by less than the Thou- 
less energy Et (related to the ergodicization time for a 
particle Tgrg by the Uncertainty Principle Et = h/Terg)- 
For mean single-particle level spacing 5, g = Et / 5- 

Since disorder breaks all the spatial symmetries, only 
time-reversal T and possibly Kramers degeneracy re- 
main. There are three classical symmetry classes|3, 
the gaussian orthogonal ensemble or GOE (T intact, 
no spin-orbit coupling), the unitary or GUE (T bro- 
ken), and the symplectic or GSE (T intact, with spin- 
orbit coupling). More recently, other classes have been 
identified for disordered superconductors Q and quan- 
tum dots constructed from two-dimensional semiconduc- 
tor heterostructures with spin-orbit couplingj^]. 

The question of how to incorporate electron-electron 
interactions in mesoscopic systems has enjoyed much at- 
tention. A recent proposal has been (somewhat inaccu- 
rately) dubbed the "Universal Hamiltonian" [a, 13 • Its 
central idea is to start with a general Hamiltonian 
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where a, s are single-particle eigenstates of the kinetic 
energy (including spin or Kramers degeneracy), and c, c^ 
are canonical fermion operators. For a given symme- 
try of the kinetic term, only special matrix elements of 
the interaction have a nonzero ensemble average|^. In 



the GOE, only V^%,^, V, 
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survive en- 



semble averaging. Matrix elements which do not survive 
ensemble-averaging are small (of typical size 5/ g) as are 
the sartiple-to-sample fluctuations of the terms which do 
survive[6J. Finally, in the large-g limit all small terms arc 



droppedJQl, leading to the Universal Hamiltonian 



where N is the total particle number, S is the total spin, 
and T — ^ cp^icp^^. In addition to the charging energy, 
Hjj has an exchange energy J and a superconducting 
coupling A. This last term is absent in the GUE, while 
the exchange term disappears in the GSE. In the large- 
g limit only interaction terms which are invariant under 
the symmetries of the kinetic term appear in Hu [^ Q . 

However, the correct way to determine if a coupling 
is relevant is to carry out a renormalization group (RG) 
analysisjg. If the coupling grows under RG, it must be 
kept, no matter how small it was initially, and irrespec- 
tive of whether it had an ensemble average in a partic- 
ular symmetry class. Using RG, the author and Harsh 
Mathur have shown that in ballistic quantum dots Hu 
is unstable to Landau Fermi liquid couplings for suffi- 
ciently strong couphngQ. Using a variant of the large- 
N approach (which subsumes the RG[ioI|). the author 
and R. Shankar have constructed [I3| a controlled the- 
ory of the quantum phase transition in the large-g limit. 
For finite g we idcntified|lCt lll| a weak-coupling regime 
which is controlled by Hjj , a strong-coupling regime with 
a distorted Fermi surface, and a fan-shaped many-body 
quantum critical reaime\i^ (QCR) in which the physics 
is dominated by collective critical fluctuations. For weak 
coupling, the QCR can be accessed at nonzero frequency 
or temperature (a;,T) > Eqcx UM^ ^^^ characteristic 
energy scale for the crossover in our case being [l^, [iJ 
Eqcx = tEt, where r is a dimensionless distance from 
the criticality. For weak coupling, perturbation theory 
around Hu is valid only for (w,T) <^ Eqcx, whereas 
the large- iV approach is valid for all {i-o^T) <^ Et\1^. 
Close to the transition (r ^ 1) the large- TV approach is 
superior, but is superfluous deep in the weak-coupling 
regime (r ~ 1). 

The focus of this paper is a class of problems where the 
QCR can be accessed at a tunable energy independent 



of, and much smaller than Ex, even at weak-coupling, 
and where the large-A'^ approach is essential to captur- 
ing the physics. These problems occur in systems where 
the single-particle Hamiltonian is undergoing a crossover 
from one symmetry class to another, say from the GOE 
to the GUEQ, acheived by turning on the orbital effects 
of a magnetic field. For Thoulcss number g 3> 1 of the 
original GOE, the crossover Hamiltonian is a, gxg matrix 



Hx{a) = Hgoe H —Hgue 
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where a is the crossover parameter. Properties 
of single eigenvector [l3|, and more recently multiple 
eigenvector [l^ correlations, have been computed in the 
crossover. For 2a^ = gx ^ 1, the following ensemble- 
averaged correlations hold for the eigenstates i^fiii), 
where fi ^ v label the states and i,j,k,l the original 
orthogonal labels: 
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The last term on the second line shows the extra cor- 
relations induced in the crossover [l^- The energy scale 
Ex = gxS/TT represents a window within which GUE 
correlations have spread. When Ex « Et, the crossover 
is complete. Similar expressions hold for the GOE to 
GSE crossover as well[l4l . 

It was pointed out|l4J that interaction matrix elements 
have enhanced sample-to-sample fluctuations in crossover 
ensembles, controlled by 1/gx and not 1/g. Thus a naive 
application of the symmetry principle leading to the Uni- 
versal Hamiltonian is problematic 
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Region II, Quantum Critical Regime 



Region I 
weak-coupling 



?Qcx=l"-"*IEx=rEx 



. '-" Region III 
„'-■ strong-coupling 



FIG. 1: Different regimes and crossover energy scales 

We are now ready to state our central results. Con- 
sider phase transitions driven by interactions which com- 
mute with the symmetries of the initial ensemble but 
not with those of the crossover ensemble. Then, (i) The 
energy scale Eqcx characterizing the many-body weak- 
coupling to quantum-critical crossover is proportional to 
Ex- (ii) The large- iV approach gives a nonpcrturbativc 
description of the physics valid for energies (w, T) <C Et 
(regions I and II in Fig^l, while pcrturbative treatments 



based on Hu are valid[Tl| only for (w, T) <IC Eqcx ~ Ex 
(region I), (iii) Tuning the single-particle crossover offers 
a powerful way to control and access the quantum-critical 
crossover. These results apply equally to diffusive and 
ballistic dots. Finally, RMT is much better controlled in 
the crossover regime Ex <IC -EtJI^ than when applied to 
the entire Thouless shell [lOlllll|. 

The tunability of the QCR is potentially very impor- 
tant for quantum dots fabricated from two-dimensional 
semiconductor heterostructures, since current samples 
appear to be on the weak-coupling side of the Stoner 
and other transitions (see references in refs.[3,[il|). An 
illustration relevant for these samples is provided by 
the crossover from the GOE to the new "intermediate" 
symplcctic ensemble [3 with small spin-orbit coupling 
in two-dimensional heterostructures. Here the physical 
crossover parameter a is a function of L/Lsoi where 
Lso is the spin-orbit scattering length. Keeping the low- 
est leading order in this ratio [a = CL^/L^q, where 
C is a factor of order unity) leads [3 to a Hamiltonian 
which conserves Sz, but not the total spin. The Univer- 
sal Hamiltonian for this case has recently been worked 
out by Alhassid and Rupp[l3| (AR). For the sake of sim- 
plicity, let us ignore all interactions except exchange. The 
Hamiltonian can be written as 
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The last term would be absent in the Universal 
Hamiltonian[l3|, since it does not commute with the ki- 
netic term, but physically it should be present because 
the hig h-energy processes which led to it are still present. 
AR|l7l| show that the spin operators can be written as 
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with S- = (S'-i-)T. The quantities Af^i, satisfy 
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where Ex = C^a'^irEx (similar to Eq. Q). It is 
straightforward to carry out the Hubbard-Stratanovich 
transformation, integrate out fermions, and obtain an 
effective action for the collective variables, here a^i <^y 
couphng to Sx, Sy. Deep in the crossover Ex ^ 5 the 
effective action is self-averaging flOl llH . Calling J = JJ, 
and taking g — > oo and setting J^ = for simplicity, to 
quadratic order we have the Euclidean action 



S, 
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The critical coupling is J* = 1, which is also the bulk 
critical value, exactly as in ref.[l^ (in fact, the quadratic 
effective action for the model of ref.[Ia|, which the au- 
thors did not compute, would be identical to Eq. ^)- 
The retarded Sj^S- correlation function in weak-coupling 



( J < 1) at T = is read off by making the replacement 
ioj ^ LO + iO^ in Eq. ||SJ|. The imaginary part of this 
correlator has the scaling form valid for oj ^ Et 

-i/m((0|5+H5_(-^)|0)) = --^ (9) 

where r = 1 — J and x = lj / Exf is the scaling vari- 
able. For Lo ^ Exr = Eqcx the system crosses over 
to the QCR. Eq. illustrates all our central results, 
the fact that Eqcx ot Ex, the fact that the large- A^ 
approachjl^ [ll| captures the nonperturbative physics of 
collective critical fluctuations, and the tunability of the 
crossover to the QCR. 

To illustrate the generality of these results, consider 
the GOE to GUE crossover in three-dimensional ballis- 
tic/chaotic superconducting nanoparticles of linear size 
L. The single-particle spacing is i5 ~ H^ /mkpL^ and the 
Thouless energy is Et ~ H^kF/mL, leading to a Thou- 
less number g « (kpL)^. A commonly-used model to 
understand superconductivity in such particles [l^ in the 
absence of an orbital magnetic field is the reduced BCS 
Hamiltonian 
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where i, s are orthogonal state and spin labels, A > 
is the attractive dimensionless BCS coupling valid in an 
energy shell of width 2tiujD around the Fermi energy, and 
T = X!ci,xCi,T- Much recent work has concentrated on 
ultrasmall grains[l8J of size a few nanometers, where 5 is 
comparable to the bulk BCS gap A = 2/iwd exp— 1/A. 
The orbital effects of an external magnetic field can be 
captured by the crossover Hamiltonian of Eq. Q , with 
"^/■y/ff = C{(t)/4>o), where 4> = BV^ is the flux through the 
sample, ^o — h/ e is the flux quantum, and C is a con- 
stant of order unity. This leads to Ex = ^^ ( 4- ) Et oc 

B^L^. For grains of Pb or Al smaller than lOnm and 
fields in the sub-Tesla range this crossover scale is smaller 
than (5, and thus orbital magnetic effects are of no im- 
portance for their physics. However, for grains in the 
size rang e lO-SOnm in moderate fields (see, for example, 
ref. [ia|)j the orbital effects will be relevant, while the 
Zeeman coupling can be ignored. We will assume that 
the critical temperature of the superconductor is smaller 
than the Thouless energy of the particle, true of Pb or Al 
particles in this size regime. Also, L should be smaller 
than the London penetration depth. In the presence of 
an orbital B field, the kinetic energy of Eq. H10|l will be 
replaced by a random matrix drawn from the ensemble of 
Eq. Q . Since the high-energy processes which led to an 
attractive BCS interaction are not modified by the small 
orbital field, the BCS interaction must be kept, though 
it no longer commutes with the symmetries of the kinetic 
energy. This leads to 
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To study the magnetization of the particle in the 
crossover we start with the partition function Z = 
Tr{exp—(3H) where [3 = 1/T is the inverse temperature. 
We convert the partition function into an imaginary time 
path integral and use the Hubbard-Stratanovich identity 
to decompose the interaction, leading to the imaginary 
time Lagrangian 

C^^-^-y^c^Adr-e^y^^s+aT + aT (13) 
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where a, a are the bosonic Hubbard-Stratanovich fields 
representing the BCS order parameter and c, c are Grass- 
man fields representing fermions. The fermions are inte- 
grated out, and the resulting action for cr, a is expanded 
to second order to obtain 
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where 



27rn//3, and the sums are restricted to 



|e^|,|ey| < hwD- We see that the correlations between 
different states /^, v play an important role. Deep in the 
crossover (for Ex ^ 5) we can replace jM^^p by its en- 
semble average (dominated by the last term of Eq. Q), 
just as in our previous work [l(l.lll] |. We will also hence- 
forth replace the summations over energy eigenstates by 
energy integrations with the appropriate cutoffs. The 
sums can be approximately carried out to give 
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where C" ~ 3.08 is chosen to obtain the correct transi- 
tion temperature in the absence of a magnetic field. The 
coefficient of the quadratic term for the iojn = term as 
T^O, /3^ooisi— log ^^"° , leading to a continuous 

phase transition at A* = 1 / \og{2hjj o / Ex) ■ 

On the normal side (A < A*) a, a fluctuate, and the 
quadratic action is sufhcient to explore their fluctuations. 
One can now integrate out cr, a to obtain the contribution 
of their thermal fluctuations to the free energy. One then 
obtains for the magnetization 
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where Mnonint is the contribution to the magnetization 
from noninteracting electrons (intimately connected to 
the noninteracting persistent current [l|). 

To see t he p hysics of the crossover into the QCR as 
T increases [12 in the most transparent way, we take the 
weak-coupling limit A — > 0, hujT) -^ 00 such that T^ re- 
mains invariant. Defining r = A^^ — (A*)~^, we obtain 
for the fiuctuation magnetization in the scaling region 



r= EAf^,c,,;c^,t: A'V=EV'/.(«)^.(i) (12) -^^0 ^ 
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where x = T/Ex is the scaling variable, and Mq = 
L'^ ^ {^C^ ExEt / T^ffi) sets the T-indcpcndcnt scale for 
the magnetization. At quantum criticality r = 0. For 
r > and T 3> Ex the second log in the denominator 
dominates, and the system crosses over into the QCR. In 
this case, unlike the previous one, the critical coupling 
A* depends on Ex, making the critical point tunable as 
well. Note that this result is nonpcrturbative in Ex, T, 
and A, and offers a perspective complementary to pre- 
vious results which are perturbativc in one or more of 
these parameters J2ri l2ll l2^ . The regime of validity of 
Eq. (PJ is £;t > (T, Ex) > (5, and r > 0. 

In summary, we have identified a class of mesoscopic 
problems in which access to the quantum critical regime 
(where the physics is dominated by collective critical 
fluctuations) can be obtained at energies Eqcx rnuch 
smaller than Et by tuning the crossover scale Ex be- 



tween ensembles with different single-particle symme- 
try. Perturbativc approaches in these systems are lim- 
ited to (y,T ) ^ Eqcx ~ Ex, whereas the large- iV 
approach|lOt lllj is successful in capturing their non- 
pcrturbative physics in the entire range {uJ,T) <^ Et- 
Symmetries crossovers in weak-coupling large-g dots of- 
fer theoretical control over and experimental access to 
strongly correlated physics. It would be interesting (and 
directly relevant to superconducting particles) to see the 
crossover from the strong-coupling regime to the quan- 
tum critical regime, as well as to determine the signatures 
of the weak-coupling to quantum critical crossover in the 
transport properties of semiconductor quantum dots. 
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